Canonical band theory of non-collinear magnetism 
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A canonical band theory of non-collinear magnetism is developed and applied to the close packed 
fee and bee crystal structures. Several examples of non-collinear magnetism in the periodic table 
are seen to be canonical in origin. This is a parameter free theory where the crystal and magnetic 
symmetry, and exchange splitting, uniquely determine the electronic bands. The only contribution 
to the determination of magnetic stability is the change in band energy due to hybridisation resulting 
from spin mixing, and on this basis we are able to analyse the origin of the stability of non-collinear 
magnetic structures, and the instability of the FM state towards non-collinear ordering. 

PACS numbers: 75.10.-b 



I. INTRODUCTION 

Structural trends in the periodic table are now very 
well understood, in contrast the situation for magnetic 
trends is not so clear. This can be attributed both to the 
fact that only a subset of the periodic table is magnetic 
in the three dimensional solid - the Fe group and rare 
earths - and also the much greater degree of freedom for 
magnetic ordering as compared to structural ordering. 
Indeed, within the transition metal block (with the ex- 
ception of Mn) all elements are found either as fee, hep, 
or bec lattices whereas a wealth of magnetic structures 
may be found. One has has both longitudinal and helical 
spin density waves, ferromagnetism, antiferromagnetism, 
as well as other more general non-collinear structures. 

For the Fe-group metals at ambient pressure the sit- 
uation is that near the centre of the series one typ- 
ically finds antiferromagnetic (AFM) structures whilst 
near the ends ferromagnetic (FM). In between one has 
several metals - 7-Fe (fee Fe), and bec and fee Mn - 
for which non-collinear structures can be found 1 . On 
the other hand, the late rare earth group elements are 
non-collinear, whilst in the middle of the series one has 
ferromagnetic Gd 2 -. 

Of course, in a fully relativistic picture the magneti- 
sation vector field will always, to some extent, be non- 
collinear. However, for the Fe-group and rare earths 
this leads only to intra-atomic non-collinearity and not 
to the coarse grained inter-atomic non-collinearity where 
the average moments associated with different sites point 
in different directions. 

The richness of magnetic structures means that most 
theoretical approaches focus on specific materials, and 
only a few attempts have been made to determine con- 
ditions whereby one places possible magnetic structures 
in some general scheme. A notable early attempt along 
such lines was the work of Pettifor 3 who, on a similar 
basis to Friedel's theory for the stability of crystal struc- 
tures, formulated a phase diagram in terms of the Stoner 



parameter to the band width I/W and the d-band oc- 
cupation number. In this theory the ferromagnetic and 
disordered local moment structures were considered, and 
the phase diagram revealed a simple argument as to why 
bec Fe is a good local moment system, but fee Co and fee 
Ni are not. Heine and Samson^, again using arguments 
based on generalised Stoner criteria constructed a similar 
phase diagram diagram but included also the AFM struc- 
ture. They further made actual calculations of the gener- 
alised Stoner I within the tight binding approximation, 
showing how 7-Fe was placed at the crossing point of the 
AFM and FM stability criteria, and hence was likely to 
assume a non-collinear structure. Hiraij^ on the basis of 
an approximation to Hartree-Fock theory, calculated the 
energy of FM, AFM, and helical spin spiral structues as 
a function of an intra-atomic parameter and the c?-band 
occupation number. The resulting phase diagram showed 
the appearence of a region of helical stability inbetween 
regions of FM and AFM stability. Unfortunately, the 
electronic structure in the Hartree-Fock approximation 
is known to be drastically in error compared to experi- 
ment. 

Recently, an attempt was made to re-examine the is- 
sue of non-collinear stability on the basis of first principle 
calculations 6 . Interestingly, it was found that materials 
collinear at their equlibrium moment were non-collinear 
for smaller (local) moments, induced either by pressure 
or a fixed spin moment procedure. On this basis, it 
was speculated both that any magnetic material could 
be made non-collinear and also that the primary factor 
governing the instability of the FM state towards a non- 
collinear state was the hybridisation of crossing spin-up 
and spin-down bands at the Fermi level. However, in 
contrast to previous approaches no attempt was made to 
link this to the ci-band occupation. 

The purpose of the present article is to unify the spirit 
of the earlier works, where the emphasis was upon find- 
ing general magnetic phase diagrams, with the approach 
in Ref. @| where the importance of spin hybridisation for 
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specific magnetic structures was emphased. As we shall 
show, this can be accomplished on the basis of Ander- 
sen's canonical (i-band theory^, suitably generalised to 
deal with non-collinear magnetism. It has been shown 
in the past that canonical <i-band theory gives a good 
qualitative account of crystal stability in the transition 
metal block and rare earths^, even including the im- 
pact of collinear magnetism upon crystal stabilit y 10 ' 11 . 
Until now, however, it has not been used to investigate 
non-collinear magnetic stability. 

The remainder of this article is structured as follows. 
In Section II we describe canonical band theory and its 
generalisation to deal with non-collinear magnetism, fol- 
lowed in Section III by details pertinent to its implemen- 
tation. In Section IV we present non-collinear magnetic 
phase diagrams for the close packed lattices, and in Sec- 
tion V discuss the reasons for non-collinear stability. Fi- 
nally, in Section VI we conclude. 

II. NON-COLLINEAR CANONICAL BAND 
THEORY 

Within the Linear Muffin Tin Orbital electronic 
structure method in the atomic sphere approximation 
(LMTO-ASA) the electronic bands are given by the sec- 
ular equation 

[PL>L(E)8 Vm , lm - S LIL {k)]a lk = (1) 

where L = (Im) are angular momentum indices, k a re- 
ciprocal lattice vector, and i a band index. This equa- 
tion has the remarkable feature that it is split into a part 
which depends solely on the one-electron potential for 
the given atomic species, Pl'l(E), known as the poten- 
tial function, and a part which depends only on the crys- 
tal symmetry, the structure matrix Sx'L(k). Neglecting 
the off diagonal VI blocks of the structure matrix, which 
amounts to neglecting all hybridisation effects between 
different angular momentum chanels, leads to a set of 
pure / canonical bands given by Pu = Su(k) unique for 
each crystal structure. In cases where the physics of the 
(i-band is expected to play a leading role one can further 
neglect all but the I' = I = 2 block and by approximat- 
ing the potential function by P = (E — C)/(/j,S 2 ) hnd a 
simple eigenvalue problem for the canonical d-bands 

[(e ik -C)/(^ 2 )l-S(k)]a lk -0 (2) 

where 1 < i < 5. One should note that since C is the 
d-band centre and l/(/zS' 2 ) simply sets the energy scale, 
this is a parameter free theory. 

To develop a non-collinear magnetic generalisation one 
first must specify a particular class of magnetic struc- 
tures. It is convenient to consider helical spin spiral 
structures, since these contain both the FM and AFM 
structures as limits. In a helical spin spiral structure the 
magnetisation is given by 



m(R„) = (sin 9 cos(q.R n ), sin 9 sin(q.R„), cos 9) (3) 

with R„ a direct lattice vector and q a reciprocal lat- 
tice vector. Although this structure may be incommen- 
surate with the underlying lattice, it is invarient under 
a lattice translation R„ and a rotation about the spiral 
axis by q.R„. This leads to a generalised Bloch theorem 
and an LMTO-ASA secular equation. By deploying the 
canonical procedure outlined above one finds an eigen- 
value problem for the canonical bands of the spin spiral 
given by 

[(e,k. - C)/(vS 2 )l - (S(k) + A)]a ik(T = (4) 

where 

cn s _ / S+ + S~ cos0 S~sin6> \ 
SilcJ -^ S-sin0 S+-S-cos9j W 

-(TJ /2 ) M 

. Eq. [4] determines a set of bands unique for each sym- 
metry of the spin space group. 

To complete the theory one needs to determine the 
canonical moment m projected onto the local quantisa- 
tion axis, for a given exchange splitting. This may be 
extracted by evaluating the expectation value of u z over 
the canonical wave function which leads to the natural 
result 

occ 
ik 

and where the sum is over all occupied states. To- 
gether with the (i-band occupation number n and the 
one-electron energy e, given respectively by 

occ 

n = ^Ifliktl 2 

ikcr 
occ 

e = ^ |aiko-| 2 £iko- (8) 

/krr 

these form the basic quantities of the theory. However, in 
order to determine the stability of competing magnetic 
structures for a given n and m what one needs is not 
the one-electron energy but the kinetic energy. This is 
due to the fact that in a canonical theory there is no in- 
trinsic intra-atomic force causing the magnetisation, but 
rather it is imposed by setting the exchange splitting to 
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some finite value which thus plays the role of an exter- 
nal field. The contribution of the energy of this external 
field is included in the one-electron energy and so must 
be subtracted, which leaves the kinetic energy of the one- 
electron system 



T = e- I A(m)dm (9) 
Jo 

as the quantity to be used for comparing the stability of 
different magnetic structures. 



III. COMPUTATIONAL DETAILS 

In order to solve Eqs. 0] to [8] we use a Monkhurst-Pack 
mesh with 18900, 109800, 216000, and 37820 k-points in 
the irreducable Brillouin zone for the FM structure, and 
spirals along FX, XW and WL, and LG respectively. For 
the bcc irreducable Brillouin zone we use 5200, 28830, 
and 109800 k-points for the FM structure and spirals 
along TH, and HP respectively. For the k-integration in 
Eqs. [7] and [8] we use an eigenvalue broadening of 4.5mRy. 
For the potential parameters which, for the canonical ri- 
band theory, simply set the origin and scale of energy we 
set C = and take l/(fj,S 2 ) for non-magnetic fee Fe. 

One should note that the natural variables of the the- 
ory outlined above are the Fermi energy Ep and the ex- 
change splitting A. However, the physical variables for 
determining magnetic stability are the ci-band occupa- 
tion number n and the spiral local moment to. Thus we 
need to solve n'(Ep, A, q) = n and m'(Ep, A,q) = m 
to find the values of Ep and A that give the required n 
and m. This can easily be done by minimising a function 
defined as 

/ = [n'(E F ,A, q) - nf + [m(E F ,A. q) - to] 2 . (10) 

In order to facilitate this it is convenient to pre-calculate 
the functions n', to', and also the corresponding one- 
electron energy function e'(Ep, A, q) over the full range 
of Ep and A and interpolate using cubic splines. In this 
way one can very rapidly minimise Eq. 1101 whereas if n! 
and to' are repeatedly calculated by solving Eq. [4] over 
the Brillouin zone this would be a very slow procedure. 



IV. CANONICAL PHASE DIAGRAMS 

In Figs. [T]and Owe show the calculated phase diagrams 
when one considers spin spirals along the TXVF symme- 
try lines of the fee Brillouin zone and the VHP lines of 
the bcc Brillouin zone. In order to construct these phase 
diagrams we have used planer spin spin spirals, 8 = tt/2 
in Eq. [31 since then AFM structures are limits of spin 
spirals. In fact, for single basis crystal structures like fee 



FIG. 1: Canonical phase diagram for FX IT spin spirals in the 
fee structure. 
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FIG. 2: Canonical phase diagram for VHP spin spirals in the 
bcc structure. 
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and bcc, all points of the BZ boundary correspond to 
AFM structures. 

One can notice immediately that, as expected, near 
half filling of the d band one finds AFM structures while 
for the nearly fully occupied or empty d-band one finds 
ferromagnetism. Unfortunately, the correspondence with 
the d electron numbers of the transition metals is only 
approximate, since the region of FM is too small in both 
the fee and bcc phase diagrams. For instance, in the fee 
phase diagram FM is the stable solution for n > 8.8, 
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whereas since fee Co is FM the boundary should be n « 
8.0. 

In Ref. la it was noted that fee Co developed a spin 
spiral structure along the TX direction if the moment 
was reduced from its equlibrium value to 0.8 fiB- I n 
fact, as one can see from Fig. [T]this is what happens in 
the canonical picture as well: there is a transition from 
FM stability to a spin spiral in the TX direction as the 
moment is reduced for n ss 8.8. Interestingly reducing 
the value of n continuously increases this splitting about 
the r point so that it becomes a spiral with q = [0, 0, 0.5], 
which is exactly the spin spiral structure observed for 7- 
Fe in calculations. In fact, in the fee phase diagram it can 
be seen that there is quite a large region where this spiral 
is stabilised. Furthermore, by lowering the moment for 
n « 6.5 this spiral closes about the X point until one 
finds nearly degenerate spin spirals along the TX and 
XW symmetry lines, again as is seen in self-consistent 
calculations. Thus one can see that the formation of 
small moment spin spirals for larger n moves continuously 
to the intermediate moment spin spiral of 7-Fe before 
finally the X AFM structure becomes stable near half 
filling. This is an appealing picture, since it unites the 
results of Ref. into a general magnetic phase diagram 
for the fee structure. In Fig. [3] are shown the spin spiral 
spectrum's corresponding to the various magnetic phases 
discussed. 

For the bee structure, shown in Fig.[2]one can note that 
the region of stability for the FM structure is much larger, 
which makes sense since bee Fe is ferromagnetic whereas 
7-Fe is non-collinear. However, as for the fee lattice the 
region of FM stability is too small. Interestingly, at n w 
9.0 one finds an incursion of non-collinear structures for 
small moments. Again, this is what is seen in Ref. 0, 
where bee Ni becomes non-collinear for small moments. 

Finally, one can note the approximate symmetry of 
both the fee and bee magnetic phase diagrams about 
half filling. This is interesting since, as was already men- 
tioned, the late rare earths all assume a non-collinear 
structure. Since for these metals the / electrons are lo- 
calised and simply play the role of polarising the itincr- 
ent spd electrons one might expect that again the cru- 
cial parameters for describing their magnetism would be 
the d-band electron number and moment m. In fact, 
this is exactly what happens for the structural energies 
where for 1.5 < n < 2.0 one finds the sequence of rare 
earth structures^. The symmetry seen in Figs. [1] and 
[2] is suggestive of the fact that there is a connection be- 
tween the late rare earths and the late Fe-group elements 
with artifically lowered moments. However, in contrast 
to the canonical crystal energy differences, the magnetic 
canonical theory is not sufficiently quantitative to give 
the actual non-collinear structures of the rare earths, and 
therefore we have not shown the hep phase diagram. 

For a more quantitative theory, the effect of angular 
momentum hybridisation must be included, but in this 
case it is no longer a parameter free theory, since band 
centres, band distortions, and band widths must be in- 



FIG. 3: Canonical spin spiral spectrum's for the symmetry 
path rXWLr in the fee Brillouin zone. Shown are spectrums 
for (a) n = 9.3, m = 0.5, (b) n = 8.8, m = 0.4, (c) n = 7.5, 
m — 1.2, and (d) n = 0.5, m = 3.3. Shown in the lowest panel 
is (e) the region around the X point for n — 6.5, m = 0.8. The 
top scale refers to graphs (a)-(d), and the bottom to graph 
(e). 
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FIG. 4: Canonical Density of States for the (a) ferromagnetic, 
(b) q = [0.0, 0.0, 0.5] spin spiral, and (c) antiferromagnetic X- 
point structure. 
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troduced for the I = 0, 1, 2 chanels. 

V. ORIGIN OF STABILITY OF 
NON-COLLINEAR STRUCTURES 

It is interesting to first reconsider the arguments for 
the stability of the collinear AFM and FM structures. 
The reason for the existance of AFM near half filling and 
FM for the nearly full band is well known and has been 
described by a number of authors (see Ref. 0] and refer- 
ences therein). The hybridisation between the different 
sublattices of the AFM structure, equivalent to the hy- 
bridisation from spin mixing when the AFM structure is 
considered as the limit of a planer spin spiral, opens a 
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large pseudo-gap in the density of states. It is this which 
stabilises the AFM structure at half filling. The argu- 
ment is thus couched in terms of significant movements 
of spectral weight at all energies, and not to changes at 
only the Fermi level. In Fig. [4] are shown the density of 
states (DOS) for the FM, X point AFM, and spin spiral 
q = [0.0,0.0,0.5] structures for the fee lattice. As can 
be seen the structure of the spiral DOS is somewhere be- 
tween that for the FM and AFM structures. Thus the 
occurance of non-collinear stability between the FM and 
AFM regions can be understood by the same reasoning 
that explains the placement of the FM and AFM regions 
themselves. 

It is natural that the determination of the lowest en- 
ergy magnetic structure will, in general, come about as 
the result of movements of spectral weight at all energies 
since generally spin hybridisation will be very significant. 
On the other hand, it is not immediately clear that this 
will be the case when one considers the instability of the 
FM state towards non-collinear structures. In Ref. it 
was speculated that such instabilities are driven by the 
hybridisation of crossing spin-up and spin-down bands 
at the Fermi level. This is, of course, the classic picture 
when one has the situation of Fermi surface nesting, as 
happens in the late rare earths. There it is the Fermi sur- 
face of the non-magnetic state which shows nesting and 
this is relevent to the final magnetic state since the itin- 
erant spd moment is very small (0.2-0.6 /is). However, 
the argument presented in Ref. ly is that the number of 
spin-up and spin-down band crossings at the Fermi en- 
ergy may more generally hold the key to understanding 
the instability of the ferromagnetic state towards a non- 
collinear state. 

The stability or otherwise of the FM state is given 
by the eigenvalues of the Hessian of the FM point in 
a general energy space of magnetic structures. For the 
magnetic structures considered in this work, this could 
be determined either as a 9 — > limit for a finite q or 
alternatively asaq^ [0, 0, 0] limit for finite 9. To inves- 
tigate the origin of the instabilty of the FM state towards 
non-collinearity we choose to investigate the energy dif- 
ference of the FM structure with the q = [0.0,0.0,0.1] 
spin spiral. This will be positive if the FM structure is 
stable, and negative otherwise. If the arguments of Ref. 6 
are relevent then this energy difference should be domi- 
nated by the contribution from a window near the Fermi 
energy. 

In Fig. [5] is plotted the energy difference of the FM 
and q = [0.0,0.0,0.1] structures for the fee lattice as 
a function of band filling for a fixed exchange splitting 
of A = 0.3. This leads to a path in the (n,m) phase 
diagram with m almost at the saturation value for each 
n (at half filling m = 4.2). Also shown in Fig. [5] is a 
breakdown of the energy into that arising from a small 
window near the Fermi energy, with AE = 0.015, and 
that from the remaining spectrum. One can clearly see 
that it is this energy and not that arising from the Fermi 
energy window which follows the total energy difference. 



FIG. 5: Energy difference of FM and q = [0.0,0.0,0.1] spin 
spiral structures. Full line is total energy different whilst cap- 
tion indicates components. 
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In fact, while the indication of FM stability from the 
total energy difference is in agreement with what one 
would expect from the phase diagram in Fig.[TJ the Fermi 
window energy difference shows no correspondence to it. 
One can also notice that the energy difference has exactly 
two crossing points. This is in line with the theorem of 
Heine and Samson 12 stating that any property of a tight 
binding band must have at least two crossing points as 
a function of band filling. A similar picture is found for 
other settings of the exchange splitting. 

It would appear then that there is no special signifi- 
cance attached to spin hybridisation at the Fermi level 
in considering the stability of the FM state. This makes 
sense when one considers that, for a given band filling, 
if the hybridisation arising from solely around the Fermi 
level dominated, then the corresponding change for some 
lower energies must integrate to zero. This could only be 
the case if the change in band structure from hybridisa- 
tion is exactly symmetric. This cannot, of course, be the 
case for all energies i.e. for all c?-band electron numbers. 
Thus, in general, it is hybridisation changes for all ener- 
gies that determine even the stability criteria of the FM 
state towards non-collinearity. 



VI. CONCLUSIONS 

To conclude a generalisation of Andersen's canonical 
band theory has been used to derive magnetic phase di- 
agrams for the close packed bec and fee structures. One 
finds that a number of the non-collinear structures ob- 
served in self-consistent calculations are found, but that 
the phase diagrams are in error for moments near satu- 
ration and the regions of ferromagnetic stability are too 
small. This is most likely due to the missing hybridisation 
between different angular momentum chanels. A similar 
situation arises in the case of crystal stability where er- 
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rors found in the d-band canonical theory are removed 
by inclusion of s and p states^. 

An analysis of the canonical magnetisation energetics 
in terms of contributions from different energy windows, 
shows that, generally, it is hybridisation from the whole 
spectrum which acts to determine the stability of the fer- 
romagnetic state with regard to non-colliear magnetism. 
Thus the picture of band crossing at the Fermi level em- 
phased in Ref. H does not appear to be relevant for the 
determination of the stability of the ferromagnetic state. 
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